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ABSTRACT 

The phase transition kinetics of Ising gauge models are investigated. Despite the 
absence of a local order parameter, relevant topological excitations that control 
the ordering kinetics can be identified. Dynamical scaling holds in the approach 
to equilibrium, and the growth of typical length scale is characteristic of a new 
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universality class with L{t) ~ (t/lnt) ' . We suggest that the asymptotic kinetics 
of the 2D Ising gauge model is dual to that of the 2D annihilating random walks, 
a process also known as the diffusion- reaction A + A ^ inert. 



PACS numbers: 05.50.+q, 05.40.+j, 05.70.Ln 



In the approach to equihbrium, systems undergoing phase transitions typi- 
caUy exhibit certain degrees of statistical self- similarity, with an example being 
the occurrence of power-law domain growth and dynamical scaling of the corre- 
lation functions!^' These scaling behaviors are believed to be quite generic, and 
can be roughly categorized into a few universality classes. The dominant factors 
affecting the classification, among others, include spatial dimensionality, conser- 
vation property, and the order parameter symmetry. Most studies in this area 
have focused on systems involving spontaneous symmetry breaking. In such sys- 
tems, it has been recognized that topological defects, i.e., singularities in the local 
order parameter, play an important role in controlling the long-time dynamics. 
In fact, all major features of the growth kinetics can be interpreted in terms of 
the behavior of these topological defects 

In this paper, we investigate the ordering dynamics in the kinetic Ising gauge 
models, which, even though exhibiting genuine phase transitions, do not have a 
local order parameter. These systems are unusual since they are invariant under 
a local gauge symmetry which can not be spontaneously broken!^' This work is 
partly motivated by the recent observation'"'' that systems involving gauge fields 
show interesting new features in dynamics. We shall illustrate that, although 
the absence of a local order parameter precludes the characterization of topolog- 
ical defects using the usual homotopy analysis, it is still helpful to find the low 
energy topological excitations that are the disordering agents of the system. In 
the Ising gauge models, they are identified as frustrated plaquettes or frustrated 
plaquette-strings. These objects, being gauge invariant, dominate the ordering 
kinetics asymptotically. While dynamical scaling is still found to hold during 
phase transition, the coarsening of typical domain size follows a unique growth 
law L(t) ~ (t/lnt)^^'^ in both two- and three-dimensions. We suggest, in partic- 
ular, a duality between the asymptotic kinetics of the 2D Ising gauge model and 
that of 2D annihilating random walks. 

Consider an infinite d-dimensional hypercubic lattice with Ising spins u = ±1 
defined on each link. The spins couple ferromagnetically to each other, and are 
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all coupled to a heat bath at temperature T. Let the pair (r, a) label a link with 
r being the coordinate of a site and a denoting a unit lattice vector. Equilibrium 
properties of the system are uniquely determined by assuming the Hamiltonian 
to take the Ising form: 

r,a,/3 

where the four-spin products are formed around the plaquettes of the lattice. 

This Ising gauge model differs from the usual Ising model in that it has a 
local symmetry while the latter has a global symmetry. Indeed, the Hamiltonian 
(1) is invariant under the gauge transformation ar,a 9rO'r,a9r+a with = 
±1, corresponding to the local operation of flipping all the spins connected to 
a lattice site. The above model belongs to a more general class of lattice gauge 
models first proposed by Wegner!*' Its equilibrium properties are well-known!"' In 
two dimensions, exact duality exists between the partition function of the Ising 
gauge model and that of the uncoupled ID Ising models. Thus, the 2D Ising 
gauge model has a phase transition strictly at T = 0. Another duality between 
the 3D Ising gauge model and the 3D Ising model proves the existence of a 
finite-temperature second order phase transition. For d > 4, numerical evidence 
suggests that the model has first order phase transitions at T 0. 

Dynamics are incorporated into the model by assigning a transition rate for 
single spin flip ar,a — ^ — Cr,a- We use a rate of the Glauber form!'^' 

W{ar,a -(Tv,a) = ^ " ^r,a tanh E{{a}r,a)] (2) 

where {o'}r,a denotes the set of all spins except crr,a ^ind 

E{{a}r,a) = {J/kBT) ^ Crr+a,l3Crr+p,aCrr,p (3) 

with the P summation restricted to unit vectors distinct from a. Rate (2) satisfles 
detailed balance. 
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We first study the coarsening dynamics of the Ising gauge model in two- 
dimensions. Since the phase transition occurs strictly at zero temperature, we 
consider the dynamics of the system following a rapid quench from the disordered 
state at T = oo to T = 0. To proceed, we must identify the proper disordering 
agents in the phase transition. It is easy to see that the gauge model has an 
infinitely degenerate ground state at T = 0. In fact, any spin configuration 
that is a pure gauge a"r,a = g^gv+a is a ground state. In the ground state, any 
four-spin product around a plaquette has to be unity: uuuu = 1. Apparently, 
the lowest energy topological excitation of the ground state corresponds to the 
spin configuration with a "frustrated" plaquette, around which the product of 
spins equals aaaa = — 1. The frustrated plaquettes, which are themselves gauge 
invariant objects, play the roles of topological defects in governing the asymptotic 
ordering kinetics. The dynamics of one frustrated plaquette is easy to visualize 
once we realize that (2) takes a simpler form at T = 0, 



Imagine one frustrated plaquette and its neighborhood. By flipping any one of 
the four spins around the plaquette according to (4), the plaquette can either 
move to a neighboring position or annihilate with another frustrated plaquette 
near by. Plus, at T = 0, no new frustrated plaquettes will be created during 
the phase transition. Thus the ordering dynamics is a continuous process of the 
motion and annihilation of the frustrated plaquettes remnant from the initial 
state. A characteristic length scale of the system, can consequently be 

introduced as the average distance between frustrated plaquettes. 

The proceeding discussions hints at a possible connection of the current prob- 
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lem to the problem of 2D annihilating random walks (ARW). Indeed, we suggest 
that the asymptotics of the 2D kinetic Ising gauge model are identical to the lat- 
ter problem, once we regard the plaquettes as diffusing, "point" -like particles. 
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Our claim is supported by numerical simulations of the gauge model detailed 
below. It is worth pointing out that this relation between the two problems is 
an extension of the well-known duality between the zero temperature ID kinetic 
Ising model and the ID annihilating random walks!^' A n-dimensional ARW, also 
known in chemical kinetics as the diffusion- limited reaction^"' A + A ^ inert, 
is defined as the process of identical particles performing simple random walks 
and annihilating upon encounter. The ARW, as one of the simplest interacting 
particle systems,'"' has been intensively studied. One quantity of interest, the 
asymptotic mean density p{t) of remaining particles in n-dimensions, has been 
rigorously shown to behave as 

pit) ~ — as t — > oo n = 2, 

' t ' (5) 

~ as t — > oo n > 3. 

The presence of a logarithmic term forn = 2 reflects the relevance of spatial fluc- 
tuations in the walker-density'^"' and is intimately related to the weak recurrence 

[12] 

property of 2D simple random walks. 

Returning to the gauge model, we introduce a nonlocal, gauge-invariant, 
operator measuring the density of frustrated plaquettes: 

Atx = ^ (1 - craaa) (6) 

where the four-spin product is around a plaquette centered at position x. The 
average density of frustrated plaquettes is then given by p{t) = {px}- We have 
carried out numerical simulations of the 2D gauge model on square lattices of 
size 1000 X 1000. In Fig.l, we show the gauge invariant density of frustrated 
plaquettes, obtained from averages over 10 sets of random initial conditions, as 
a function of time. To extract the asymptotics of p{t), it is beneficial to plot the 
quantity tp/logt against 1/logt, as shown in the inset. The reason for plotting 
the data in such a way is to exclude the possibility of a pure power law decay 
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without the logarithmic term. Clearly p(t) ~ 1/t is excluded since this would 
lead to a curve approaching the origin at smaller 1/logt. Similarly asymptotic 
behaviors of the type p{t) ~ l/t*^ with a ^ 1 are also inconsistent with the data 
shown in the inset. Therefore we conclude from Fig.l that tp/logt — > const > 
as t — > oo, which implies that 

p{t) ~ — as t — > OO, (7) 

asymptotically, just as in the 2D ARW. As a result, the characteristic length 
scale L{t) ~ p-V2 ^ {t/lnt)^^. 

Further insights about the ordering process can be gained by studying the 
proper correlation functions. Note that, just like in equilibrium, averages of 
all non-gauge-invariant quantities vanish. We concentrate on the plaquette- 
plaquette correlation function defined as C(x — x',t) = (//xA*x')- Fig. 2 is the 
numerically calculated correlation function. It is found that the correlation func- 
tion satisfies dynamical scaling of the form 

C(x-x',t)=p2(t)F(^^^^ (8) 

where F is a universal function. For comparison, we have also computed the 
two-particle correlation function of the ARW, which, if superimposed on Fig. 2, 
would be indistinguishable from that of the gauge model. These calculations lend 
further support to our claim that the 2D gauge model is in the same universality 
class as the ARW. No analytic determinations of F{x) is yet available. However, 
we may determine its short-distance behavior rather easily. Starting from the 
general master equation for the probability density of spins and using the rate 
(4), we derive the following equation for the plaquette moments: 

d 1 

a a 

Equations for higher moments in general show a more complex hierarchy. After 
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translational average using space homogeneity, we have 



Using (7), we obtain 



|p(t) = -4C(l,t). (10) 



C(l,t)-^ ast^oo. (11) 



Assuming dynamical scaling of the form (8), Eq.(ll) gives the short distance 
behavior of the scaling function 

F{x) ~ --^ for x^O. (12) 
ma; 

At long distances, we have trivially F(a;) = 1 as a; — > oo. 

Now we proceed to study the coarsening dynamics in three dimensions. It is 
still obvious that frustrated plaquettes disorder the ground state. However, there 
is an additional feature pertaining to the 3D system. Consider the six plaquettes 
enclosing a unit cube. The product of frustrations of these plaquettes is unity 

6 

l[iaaaa)i = 1, (13) 
1=1 

since each spin is multiplied twice. Therefore the number of frustrated plaquettes 
enclosing each unit cube has to be even. This constraint introduces a nice conti- 
nuity property that the frustrated plaquettes have to connect to form string-like 
structures. Here for each frustrated plaquette we define a unit string segment 
passing through the center of the plaquette. These "strings", either forming 
a loop or extending to infinity, are the topological excitations determining the 
asymptotics of the 3D ordering kinetics. 

Let p{t) denote the total length of the string network and define the charac- 
teristic length scale of the system L{t) as the average inter-string spacing such 
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that L{t) = 1/ \/ p(t). In Fig. 1(b), we show the behavior of pit) obtained from 
numerical simulations on a 64^ lattice, averaged over 60 initial conditions. Again, 
we find that 

pit) ~ — (14) 

which implies that the typical length scale L{t) ~ (t/lnt)^/^. The presence of 
a logarithmic term in 3D coarsening is somewhat unexpected. It suggests that 
the ordering dynamics in this case is not driven purely by the curvature of the 
plaquette-strings, but that spatial fluctuations in the string density are also rel- 
evant. We have not studied higher-dimensional orderings, but it is reasonable to 
believe that this result will remain valid for all (i > 2. 

The model we considered belongs to a more general category of kinetics Ising 
gauge models, M^^^ (cJ > n), using Wegner's notation!'^' In model M^^^, the Ising 
spins are defined on all the unit (n — l)-cubes of an infinite cZ-dimensional hyper- 
cubic lattice, and the Hamiltonian (1) is replaced by a summation of the 2n-spin 
products around unit n-cubes. Dynamics are again incorporated using Glauber's 
rates. The model we considered in this paper corresponds to n = 2. Now, extend- 
ing our previous arguments to these more general models, wc reach the general 
conclusion that the asymptotic kinetics of the zero temperature M„„ model is 
dual to that of the n-dimensional annihilating random walks, for all n > 1. Con- 
sequently, the density of frustrated n-cubes is, according to (5), p{t) ~ t~^ for 
n > 3. Defining the characteristic length scale of the system accordingly, as 
L{t) ~ p~^^^i we conclude that the system coarsens as L{t) ~ t^^'^. This is true 
for all d > n > 3. 
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FIGURE CAPTIONS 



1) Average density of frustrated plaquettes p{t) in two- and three-dimensions. 
Curves in the inset are plotted in the way, tp/ logt vs. 1/ logt, to iUustrate 
the asymptotic behavior p ~ logt/t as t ^ oo. 

2) Dynamical scaling of the plaquette-plaquette correlation function in two 
dimensions. Data are averaged over 64 initial conditions on 512 x 512 
lattices. 



10 



